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FIGURES 


The  magnet os pheric  locations  of  the  three  satellites  is 
shown  at  the  time  of  the  onset  of  the  drift-echo  event. 

The  drift-echo  energy  spectrum  is  shown  as  measured 
at  1976-059  and  SCATHA . 

The  temporal  history  at  the  onset  of  the  drift-echo 
event  of  the  count  rates  in  nine  channels  is  shown  as 
observed  at  SCATHA..... . . . 

A  plot  of  the  time  delay  between  the  two  helium 

channels  and  thrae  proton  channels  is  plotted 

versus  the  reciprocal  of  the  proton  energy . 


I.  INTRODUCTION 


The  aeen  charge  states  of  heavy  ions  trapped  in  the  earth's  magneto¬ 
sphere,  in  particular  those  of  the  abundant  heavy  ions  He,  C,  and  0,  are 
important  observational  parameters  predicted  by  various  theories  of  the  origin 
and  evolution  of  the  aagnetospherlc  plasma  (e.g.,  Spjeldvlk  and  Fritz,  1978a, 
b).  Up  to  the  present  tlae  satellite  instrumentation  capable  of  determining 
the  charge  state  of  heavy  ions  with  energies  above  ~  50  keV  has  not  been 
flown.  Consequently  there  is  motivation  to  develop  Indirect  methods  of  deter¬ 
mining  the  charge  state  of  energetic  aagnetospherlc  ions,  even  if  such  methods 
do  not  have  universal  applicability.  One  such  method  is  to  measure  the  drift 
speed  of  an  Ion  In  the  magnetic  field  of  the  earth. 

The  angular  velocity  associated  with  the  azimuthal  drift  of  an  ener¬ 
getic  Ion  Is  given  by  the  relationship  (cf.,  Schulz  and  Lanzerotti,  1974) 

Qj  «  f(ao)LE/q  (1) 

where  K  is  the  kinetic  energy  of  the  ion,  q  is  the  charge  state,  and  <xq  is  the 
equatorial  pitch  angle.  Here  it  is  assumed  that  the  energy  of  the  ion  is 
sufficiently  high  that  the  effect  of  aagnetospherlc  electric  fields  on  drift 
is  not  significant.  The  important  feature  to  note  is  that  the  drift  velocity 
of  an  ion  at  a  given  L  value  and  pitch  angle  depends  upon  the  ratio  of  its 
total  energy  to  the  effective  charge  of  the  ion.  Thus  a  measurement  of  the 
drift  speed  and  the  energy  of  an  ion  determines  the  charge  state. 


Ic  is  necessary  to  tag  an  ion  in  some  way  in  order  to  be  able  to  measure 
its  drift  speed.  The  present  experimental  procedure  uses  observations  of 
certain  transient  enhancements  in  the  ion  fluxes  which  are  generated  when  a 
group  of  ions  is  Injected  and/or  preferentially  accelerated  at  some  localized 
point  in  the  magnetosphere  (for  example  by  a  substorm  in  the  midnight  sec¬ 
tor)  .  Ions  drift  westward  around  the  earth  and  can  be  seen  at  longitudes  well 
away  from  the  Injection  point  (Lanzerotti  et  al.,  1971;  Bellan  et  al., 
1978).  Since  the  enhancement  is  localized  in  longitude,  a  transient  pulse  of 
particles  will  be  seen.  Frequently  the  ions  return  several  times  to  the  loca¬ 
tion  of  an  observing  satellite,  hence  the  term  "drift  echoes".  Because  the 
drift  speed  of  an  ion  is  energy  dependent,  the  pulse  of  Increased  flux  is  seen 
in  different  energy  channels  of  a  satellite  detector  at  different  times. 

The  drift  speed  of  an  ion  depends  upon  the  actual  configuration  of  the 
earth's  magnetic  field  at  the  time  of  the  observations.  Protons  are  used  to 
determine  the  relationship  between  arrival  time  and  energy /charge  in  the 
actual  magnetic  field  configuration  of  the  earth  at  the  time  of  the  drift-echo 
event.  Protons  can  have  only  one  charge  state,  and  no  charge  ambiguity  can 
result.  When  the  drift  rates  of  the  heavier  ions  are  measured  at  known  ener¬ 
gies,  a  comparison  of  their  drift  rates  with  that  of  protons  gives  a  measure 
of  the  ionic  charge. 

Drift-echo  peaks  frequently,  are  observed  to  have  a  complex  shape  and  do 
not  exhibit  a  simple  rise  and  fall  (Belian  et  al.,  1982),  One  sees  several 
"peaks"  very  close  together.  These  "structured  peaks”  could  result  from  sev¬ 
eral  causes,  including  an  injection  multiple  in  time  or  longitude;  there  being 
more  than  one  relatively  abundant  charge  state  of  a  given  ion;  multiple  ion 
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species  being  detected  in  a  single  sensor;  and  charge  exchange  causing  the 
charge  state  of  an  individual  ion  to  vary  significantly  on  the  drift  time 
scale.  Observation  of  the  evolution  of  a  structured  peak  with  longitudinal 
drift  can  help  to  distinguish  between  some  of  these  possibilities. 


II.  RESULTS 


The  results  discussed  in  this  report  are  from  data  acquired  by  Los  Alamos 
Instruments  aboard  the  synchronous  altitude  satellites  1976-059  and  1977-007 , 
end  from  Aerospace  instruments  aboard  the  SCATHA  (1979-007)  near-synchronous 
satellite.  SCATHA  is  in  an  orbit  with  an  apogee  of  »  43200  km  (R  -  7.8  Rg), 
perigee  of  ■  27500  km  (R  ■  5.3  Rg),  and  an  inclination  of  7.8°.  The  Los 
Alamos  sensors  of  interest  covered  the  proton  energy  ran  from  145  keV  to  800 
keV  (13  channels).  The  Aerospace  sensors  covered  the  er  jy  range  from  54  keV 
to  717  keV  for  protons  (4  channels),  from  396  keV  to  9  j  .V  for  helium  ions 
(2  channels),  and  from  1264  keV  to  2704  keV  for  CNO  ions  (3  channels).  The 
Loa  Alamos  instruments  are  described  in  a  paper  by  Baker  et  al.  (1979)  and  the 

Aarospace  instruments  in  a  paper  by  Blake  and  Fennell  (1981). 

♦ 

In  this  report  a  representative  drift-echo  event  will  be  examined  which 
was  first  seen  at  the  location  of  SCATHA.  at  ■  21:10  UT  on  25  February  1979. 
At  this  time  the  satellites  were  located  as  shown  in  Figure  1.  SCATHA  was 
between  1976-059  and  1977-007  in  longitude,  and  was  also 'at  L  ■  6.6.  All 
three  satellites  were  in  the  post-noon  quadrant.  This  chance  configuration  is 
Ideally  suited  for  study  of  a  drift-echo  event. 

Figure  2  presents  the  proton  spectrum  as  observed  at  1976-059  and  at 
SCATHA.  No  normalization  of  the  data  has  been  done.  The  spectrum  is  gen¬ 
erated  by  taking  the  fluxes  observed  in  each  energy  channel  at  the  time  of  the 
maximum  of  the  first  peak  in  each  channel.  -  Thus  the  spectrum  is  not  the 
energy  spectrum  as  seen  at  any  instant  of  time  at  the  observing  satellites. 
It  should  be  the  proton  energy  spectrum  at  the  injection  point  following 
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Figure  1.  The  raagnetospheric  locations  of  the  three  satellites  is  shown 
the  time  of  the  onset  of  the  drift-echo  event. 
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Figure  2.  The  drift-echo  energy  spectrum  is  shown  as  measured 
SCATHA. 
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acceleration  before  drift  dispersion  begins  if  strong,  energy 'dependent  loss 
processes  have  not  distorted  the  energy  spectrum  between  the  lnjectloo  point 
and  observation  locations.  The  fact  that  repeated  echoes  (multiples  of  360* 
of  drift)  are  seen  in  this  event  suggests  that  substantial  losses  did  not 
•occur  between  the  Injection  point  and  the  first  observations  of  the  event, 
which  Involved  a  drift  of  less  than  180*.  There  is  excellent  agreement 
between  the  two  sets  of  measurements  and,  furthermore,  the  data  from  the  third 
satellite,  1977-007,  fall  on  the  same  curve  (but  are  not  plotted  to  avoid 
cluttering  the  figure).  If  strong  particle  losses  were  occurring  a  sub¬ 
stantial  difference  in  the  energy  spectra  would  be  expected  at  the  three 
different  longitudinal  locations. 

A  cross-correlation  analysis  was  performed  between  the  data  from  the  290- 
340  keV  channel  and  the  400-480  fceV  channel  of  both  1976-059  and  1977-007. 
The  expectation  was  that  the  peak  in  the  lower-energy  channel  would  lag  the 
peak  in  the  higher-energy  channel  by  a  larger  time  Interval  at  1977-077  than 
at  1976-059;  1977-007  is  66*  further  than  1976-059  from  the  putative  injection 
region  (Bellan  et  al. ,  1978;  Baker  et  al.,  1979)  in  the  midnight  sector  (cf.. 
Figure  1).  The  results  of  this  analysis  are  as  follows: 


1976-059:  At  ■  0.85  ±  0.05  minutes 


1977-007:  At  •  2.05  ±  0.05  minutes 


(2) 


A  relative  shift  of  significantly  more  than  A  0.05  minutes  from  the 
given  At  rapidly  reduces  the  correlation  between  two  energy  channels,  and  thus 


Is  s  measure  of  Che  uncertainty  in  the  derived  value  of  At.  These  obser¬ 
vations  lead  to  a  measure  of  the  separation  rate  per  degree  of  longitude: 


(AQ)”1  -  (1.82  *  0.15)  x  10-2  minutes/deg 


A  dipolar  calculation  (Schulz  and  Lanzerotti,  1974)  of  the  expected  drift 


speeds  of  290  keV  and  400  keV  protons  yields 


1/0(290)  -  6.36  x  10“*  minutes /deg 
1/0(400)  -  4.61  x  10  2  minutes/deg 


(AO)  *  «  1.75  x  10“^  minutes /deg 


This  value  is  in  excellent  agreement  with  the  observations,  cf.,  (3),  and 
indicates  that  the  dipole  approximation  is  a  good  one  for  ordering  the  obser¬ 
vations  in  this  event. 

The  count  rates  as  a  function  of  time  for  selected  ion  channels  of  SCATHA 
are  shown  in  Figure  3.  These  data  are  for  pitch  angles  of  90*  ±  30*.  The 
drift  speed  of  an  ion  is  a  function  of  pitch  angle;  however,  a  dipolar  calcu¬ 
lation  predicts  only  a  52  difference  between  the  extremes  of  60*  and  90*.  The 
dispersion  in  arrival  time  as  a  function  of  ion  energy  can  be  seen  clearly  in 
Figure  3  although,  because  of  the  averaging  of  the  data  that  was  done  to  gen¬ 
erate  the  figure,  it  cannot  be  used  for  quantitative  timing  purposes.  Note 
that  the  peak  in  Che  363-717  keV  proton  channel  occurs  prior  to  either  of  the 
helium  peaks.  There  was  only  one  CN0  count  in  the  five  hours  preceding  the 
event,  and  it  was  in  the  lowest-energy  channel;  these  CN0  counts  can  be  asso¬ 
ciated  with  the  drift-echo  event  with  confidence. 


21:00  21:30  22:00 

UT 


Figure  3.  The  temporal  history  at  the  onset  of  the  drift-echo  event  of  the 
count  rates  In  nine  channels  Is  shown  as  observed  at  SCATHA. 


In  order  to  sake  a  quantitative  estimate  of  the  B/q  values  for  the  396- 
S48  kaV  and  548-960  keV  helium  channels,  a  cross-correlation  analysis  was 
performed  between  these  two  channels  and  the  104-189  keV,  189-363  keV,  and 
363-717  kaV  proton  channels.  The  results  are  shown  in  Figure  4,  where  the 
time  differences  between  the  arrival  of  helium  ions  (in  both  energy  channels) 
and  the  arrival  of  protons  (in  three  energy  channels)  are  plotted  against  the 
reciprocal  of  the  energy  of  the  proton  channel  in  question.  The  vertical  bars 
Indicate  the  timing  uncertainty.  A  straight  line  is  expected;  the  intersec¬ 
tion  of  a  straight-line  fit  with  zero  time  difference  gives  the  E/q  of  the 
helium  ions.  The  results  are  190  keV/q  and  270  keV/q. 

These  calculated  values  may  be  compared  with  the  actual  channel  thresh¬ 
olds.  For  q  •  1,  E/q  amounts  to  396  keV/q  and  548  keV/q,  and  for  q  ■  2,  E/q 
amounts  to  198  keV/q  and  274  keV/q,  respectively.  The  straight-line  fits  gave 
values  of  E/q  only  8  keV  and  4  keV  different  from  the  known  values  for  q  * 
2.  Clearly  the  helium  ions  observed  in  this  drift-echo  event  were  stripped. 

An  upper  limit  to  the  fraction  of  singly  charged  ions  can  be  derived  by 
assuming  all  the  counts  around  the  expected  time  of  arrival  of  He*,  prior  to 
the  arrival  of  He**,  are  due  to  singly  charged  helium.  The  fraction  is 

He*/He**  <  0.03  (6) 

A  low  heliur  ion  count  rate  was  seen  for  a  few  hours  prior  to  the  drift- echo 
event.  Indicating  a  pre-existing  flux.  The  count  rate  at  the  expected  time  of 
the  He4*  drift  echo  is  no  greater  than  the  pre-existing  flux,  and  thus  the 
dbeervations  are  consistent  with  there  being  no  He*  ions  in  the  drift-echo 


event 
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/  190  keV 
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Figure  4.  A  plot  of  the  time  delay  between  Che  two  helium  channels  and  three 
proton  channels  Is  plotted  versus  the  reciprocal  of  the  proton 
energy*  The  open  points  are  the  proton  correlations  with  the  548- 
960  keV  helium  channel;  the  solid  points  are  the  proton  correla¬ 
tions  with  the  396-548  keV  helium  channel. 


A  similar  cross-correlation  analysis  Is  not  practical  for  the  CNO  chan¬ 
nels  because  of  the  paucity  of  counts,  cf..  Figure  3.  Although  a  quantitative 
analysis  would  be  difficult,  several  features  of  the  data  are  noteworthy: 

1*  The  first  CNO  count  does  not  appear  until  ~0.5  minutes  after  the 
Increase  In  count  rate  In  the  396-548  keV  helium  channel.  Since  It  was  shown 
that  q  ■  2  for  helium,  the  largest  value  of  E/q  In  the  396-548  keV  pass band  Is 
548/2  *  274  keV /charge.  This  result  Implies  E/q  for  the  CNO  ions  Is  less  than 
274  keV/charge. 

2.  The  three  CNO  energy  channels  show  the  expected  energy  dispersion; 
the  counts  appear  In  order  of  decreasing  energy. 

3.  The  CNO  “peak"  Is  broader  chan  would  be  expected  if  q  equaled  1  or  2. 

4.  The  CNO  peak  la  located  ~  3  minutes  later  than  the  predicted  position 
for  q*l  as  derived  from  the  proton  observations. 

5.  There  are  a  total  of  14  counts  In  the  CNO  "peak"  (summed  In  all  three 
channels).  In  order  to  Improve  the  statistics  the  CNO  counts  have  been  exam¬ 
ined  at  all  pitch  angles,  and  not  just  In  the  60-120°  antisolar  pitch-angle 
bln  which  Is  plotted  In  Figure  3.  There  are  60  additional  counts  for  a  total 
of  74.  Not  one  appears  prior  to  the  onset  Indicated  In  Figure  3. 

These  observations  are  consistent  with  the  CNO  charge  state  being  5  or 
more.  The  probability  that  the  charge  state  was  1  or  2  and  that  the  observed 
timing  was  a  result  of  statistics  Is  nil.  For  charge  state  3  and  especially 
charge  state  4  the  observed  timing  possibly  could  be  a  result  of  statistical 
fluctuations.  The  probability  would  be  difficult  to  assess  and  likely  would 
require  a  numerical  simulation. 


III.  DISCUSSION 


The  results  of  Chls  study  show  that  the  helium  Ions  were  stripped,  and 
that  the  CNO  Ions  were  probably  of  charge  state  5  or  higher,  and  definitely 
not  charge  state  1  or  2.  These  results  Indicate  that  the  source  of  the  ener¬ 
getic  Ions  was  not  the  Ionosphere.  Furthermore  these  Ions  could  not  have  been 
resident  In  the  synchronous  altitude  region  for  a  long  time  before  acceler¬ 
ation.  If  they  had  been,  then  charge  exchange  would  have  transformed  the 
stripped,  or  nearly  stripped,  solar-wind  ions  (Spjeldvik  and  Fritz,  1978a; 
1978b)  to  a  lower  charge  state  than  Is  observed.  A  model  in  which  plasma 
sheet  Ions,  originally  from  the  solar  wind,  are  brought  In  from  the  tall  and 
accelerated  would  fit  the  observations  (Baker  et  al.,  1982). 

Fritz  and  Wllken  (1976)  have  shown  that  the  energetic  heavy-ion  fluxes  In 
the  synchronous  altitude  region  are  highly  time  variable  and  that,  above  a  few 
hundred  keV,  CNO  Ions  are  the  most  abundant.  Their  results,  and  the  present 
ones,  suggest  that  the  most-energetic  ions  in  the  synchronous  altitude  region 
result  from  Injections  of  plasma  sheet  Ions  accelerated  by  electric  fields. 
If  the  plasma  had  a  steap  energy  spectrum  prior  to  acceleration,  the  several¬ 
fold  Increase  In  the  CNO  energy  relative  to  that  of  protons,  because  of  their 
high  charge  state,  could  make  than  most  abundant  In  the  energized  plasma 
population. 
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LABORATORY  OPERATIONS 


The  Laboratory  Operations  of  The  Aerospace  Corporation  Is  conducting  exper¬ 
imental  and  theoretical  Investigations  necessary  for  the  evaluation  and  applica¬ 
tion  of  scientific  advances  to  new  military  space  systems.  Versatility  and 
flexibility  have  been  developed  to  a  high  degree  by  the  laboratory  personnel  in 
dealing  with  the  many  problems  encountered  In  the  nation's  rapidly  developing 
space  systems.  Expertise  in  the  latest  scientific  developments  16  vital  to  the 
accomplishment  of  tasks  related  to  these  problems.  The  laboratories  that  con¬ 
tribute  to  this  research  are: 

Aerophyslca  Laboratory:  Launch  vehicle  and  reentry  aerodynamics  and  heat 
transfer,  propulsion  chemistry  and  fluid  mechanics,  structural  mechanics,  flight 
dynamics;  high-temperature  ehenaomechsnics,  gas  kinetics  and  radiation;  research 
In  environmental  chemistry  and  contamination;  cv  and  pulsed  chemical  laser 
development  including  chemical  kinetics,  spectroscopy,  optical  resonators  and 
beam  pointing,  atmospheric  propagation,  laser  effects  and  countermeasures. 

Chemistry  and  Physics  Laboratory :  Atmospheric  chemical  reactions,  atmo¬ 
spheric  optics?  light  scattering,  state- specific  chemical  reactions  and  radia¬ 
tion  transport  in  rocket  plumes,  applied  laser  spectroscopy,  laser  chemistry, 
battery  electrochemistry,  space  vacuum  and  radiation  effects  on  materials,  lu¬ 
brication  and  aurface  phenomena,  thermionic  emission,  photosensitive  materials 
and  detectors,  atomic  frequency  standards,  and  bloenvlronmental  research  and 
monitoring. 

Electronics  Research  Laboratory:  Microelectronics,  CaAs  low-noise  and 
power  devices,  semiconductor  lasers,  electromagnetic  and  optical  propagation 
phenomena,  quantum  electronics,  laser  communications,  Ildar,  and  electro-optics; 
communication  sciences,  applied  electronics,  semiconductor  crystal  and  device 
physics,  radiometric  Imaging;  mil lime ter- wave  and  microwave  technology. 

Information  Sciences  Research  Office:  Program  verification,  program  trans¬ 
lation^  performance-sensitive  system  design,  distributed  architectures  for 
spaceborne  computers,  fault-tolerant  computer  systems,  artificial  Intelligence, 
and  microelectronics  applications. 

Materials  Sciences  Laboratory:  Development  of  new  materials:  metal  matrix 
composites,  polymers,  and  new  Forms  of  carbon;  component  failure  analysis  and 
reliability;  fracture  mechanics  and  stress  corrosion;  evaluation  of  materials  in 
space  environment;  materials  performance  In  space  transportation  systems;  anal¬ 
ysis  of  systems  vulnerability  end  survivability  In  enemy-induced  environments. 

Space  Sciences  Laboratory;  Atmospheric  and  Ionospheric  physics,  radiation 
from  the  atmosphere t  density  and  composition  of  the  upper  atmosphere,  aurorae 
and  airflow;  magnetospherlc  physics,  cosmic  rays,  generation  and  propagation  of 
plasma  waves  in  the  magnetosphere;  solar  physics.  Infrared  eseronomy;  the 
effects  of  nuclsar  explosions,  magnatlc  storms,  and  solar  activity  on  the 
earth's  atmosphara.  Ionosphere,  and  magnetosphere;  the  effects  of  optical, 
electromagnetic,  and  partlculets  radiations  In  space  on  space  systems. 


